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ABSTRACT 

The Pan Andromeda Archeological Survey (PAndAS) CFHT Megaprime survey of the M31-M33 
system has found a star stream which extends about 120 kpc NW from the center of M31. The great 
length of the stream, and the likelihood that it does not significantly intersect the disk of M31, means 
that it is unusually well suited for a measurement of stream gaps and clumps along its length as a 
test for the predicted thousands of dark matter sub-halos. The main result of this paper is that the 
density of the stream varies between zero and about three times the mean along its length on scales 
of 2 to 20 kpc. The probability that the variations are random fluctuations in the star density is 
less than 10 -5 . As a control sample we search for density variations at precisely the same location 
in stars with metallicity higher than the stream, [Fe/H] = [0, -0.5] and find no variations above the 
expected shot noise. The lumpiness of the stream is not compatible with a low mass star stream in 
a smooth galactic potential, nor is it readily compatible with the disturbance caused by the visible 
M31 satellite galaxies. The stream's density variations appear to be consistent with the effects of a 
large population of steep mass function dark matter sub-halos, such as found in LCDM simulations, 
acting on an approximately 10 Gyr old star stream. The effects of a single set of halo substructure 
realizations are shown for illustration, reserving a statistical comparison for another study. 
Subject headings: dark matter; Local Group; galaxies: dwarf 



1. INTRODUCTION 

Galactic halos formed from LCDM initial conditions 
in n-body simulations have approximately 10% of their 
mass in orbiting sub-halos (Diemand, Kuhlen & Madau 



2007 
and 



Springel et al. 



20081). The gravitational stirring 
reatmg ot a galactic disk and halo star clusters by 
dark sub-halos is tempered by the large numbers, high 
random velocities and broad distribution of sub-halos 
throughout the overall dark halo. However, very low 
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velocity dispersion star streams in galactic halos are sen- 
sitive to the degree to which the dark matter in the halo is 
sub-structured into thousands of orbiting sub-halos. The 
sub-halos fold and chop the star-streams and gradually 
increase the velocity dispersion around the mean motion 
to about ~15% of the halo ci rcular velocity, typically 



30kms 1 , over a Hubble time (Ibata et al.|200 2; Siegal- 
[ Gaskins fe Valluri||2008[ |Carlberg||^(M| |Voon, Jolmston 
|fe Hogg||2010p . Consequently, cool star streams are sen 
sitive indicators of the presence of the predicted dark 
matter substructure. 

About half a dozen of the currently known Milky Way 
streams qualify as cool (most confidently, Pal 5, GD-1, 
Orphan, Archeron and Styx) that is, having local ve- 
locity dispersions below about lOkms -1 , or, width less 
than about 0.1 radian as seen from the center of the 
host galaxy. A dark matter sub-halo crossing such a cool 
stream will lead to visible disturbances. Unfortunately, 
the star count data often do not yet have sufficient local 
numbers to allow statistically significant measurements 
of density variat ions relative to the galactic foreground 
and background (Odcnkirch en et al.|2001[ [Ferguson et al. 
[2002] iMajewski et al.||^004| IChapman et al.| |20()8l IGrill- 
[rnair]|2009| |Odenkirchen et al.||2009f ~ 

The Pan Andromeda A rcheological Survey (PAndAS) 
( McConnachie et al. 2009[ ) in one fell swoop provides deep 
and uniform data to a distance of about 150 kpc from 
M31's center. The spectacular star stream north-west of 
M31, mor e than 100 kpc long, was first displayed in its 
entirety in Richardson et al. (2011). The great length of 
the stream, as well as being fairly distant from the dis- 
turbing effects of the main body of M31, make the stream 
an exceptionally interesting case for analysis of density 
variations. And, the star counts have sufficient signal to 
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Fig. 1.— The spatial distribution of the [Fe/H]= [-0.6, -2.4] red 
giant sta rs in the NW region of M 31. A full field version is pre- 
sented in |Richardson et ah] ( |20 1 1 [ ) The image is 10° across in the 
tangent projection co-ordinates, which are centered in the exact 
middle of this map. 



noise to allow reliable local surface density measurements 
of the stream. 

Multiple image gravitational lensing by a galaxy of a 
background quasar s is another probe of sub-halos. |Mao| 



& Schneider ( 1998 ) pointed out that the anomalous mix 
ratios relative to a smooth potential model that accu- 
rately predicts the locations of the images are the natu- 
ral consequence of substructure in galaxy halos . A sta- 
tistical analysis of available strong lens s ystems (|Dalal fc' 



Kochanek|2002| |Kochanek fc Dalal|2004[) assuming NFW 



density profile sub- halo profiles found that about 2% of 
the halo mass (with a very large uncertainty) is in sub- 
structure. L ensing is very sensitiv e to the central density 
profile which [Springel et al. ( 2008[ ) found to be close to an 
Einasto profile, which is shallower than the NFW profile. 
The theoretical situation has become somewhat unclear, 
since a projection of a sub-h alo rich n-body galaxy simu- 
lation ( Xu et al.|2009 2010) finds image distortions gen- 
erally smaller than those seen in strong lens systems. A 
better overall understanding and tighter limits re quires 
better data on more strong lens systems ( |Treu|[2010| ) . 

In this paper we analyze M31's NW star stream for 
the density variations that sub-halos are expected to in- 
duce. We first trace the mean centerline of the stream. 
The luminosity function and total luminosity of the main 
part of the stream are estimated as an indication of the 
progenitor system. The stream's width and deviations 
from the centerline are measured. The measurement of 
the density along the stream gives the data which we use 
to test whether the stream is smooth or lumpy. An illus- 
trative simulation is used to show the expected scale and 
degree of stream structure as the amount of substruc- 
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Fig. 2. — The color-magnitude diagram of stars in the M31 field. 
The stars to which thelDotter et al. 1(12008) isochrones assign metal- 
licities, under the assumption that they ire red giants at the dis- 
tance of M31, are masked out at [Fe/H]= [0, -0.5], [-0.9, -1.1] and 
[-1.5, -3], from right to left, respectively. 

ture increases from the visible dwarfs to the thousands 
of halos that LCDM predicts. 

2. A MAP OF THE M31 STELLAR HALO 

PAndAS has acquired g and i band images cover- 
ing nearly 400 square degrees of M31-M33 with the 
CFHT Megaprime system. The images have a minimum 
S/N— 10 at %ab — 24.5 mag. De tails of data reduction 
McConnachie et al. (|2010[) and Richardson et al 



arc m 



(|2011|. We adopt thepVLcConnachie et al.|(|2U05p distance 



of 785 kpc (DM of 24.47 mag J tor M31, although we note 
a statistically consistent, but somewhat smaller, Cepheid 
distance having a DM 24.32±0.12 mag or 731 kpc fvT 
lardell et al.||2007 ). The adopted distance indicates that 
the survey reaches to Mj ~ mag at high signal-to-noise. 
We use data up to and including January, 2010. 

Image extraction yields approximately 31 million stel- 
lar objects in the images. Matching the colors and mag- 
nitudes of the stars to iso chrones of varying m etallicity 
at a fixed age of 12 Gyr (Dotter et al. 2008) identifies 
M31 Red Giant Branch (RGB) stars using the proce- 



dures of McConnachie et al. (20101. Since we are seeing 
M31 through our own galaxy, it turns out that more 
than 80% of the stars do not match RGB isochrones at 
the adopted distance. These stars are not included in the 
map, although all stars above the completeness limit are 
useful for the measurement of geometric biases of the im- 
ages. The stars are i ndividually extinction corre cted with 
galactic dust maps (McConnachie et al. 2010). We are 
left with about 4.3 million nominal red giant branch stars 
to 24.5 mag which have been assigned metallicities in the 
range of [Fe/H] =[-3, 0]. Of these stars, approximately 
25% are in the NW quadrant which is our area of interest. 
The galactic stars that happen to fall into the range of 
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colors and magnitudes that are appropriate for M31 red 
giants create a slowly varying foreground which will be 
removed in the analysis of the stream density. Note that 
a shorter distance would make the inferred luminosities 
smaller and cause the isochrone match to systematically 
yield larger metallicities. To a first approximation this 
simply offsets the metallicity range we analyze and does 
not otherwise affect the analysis. 

Star streams are identified as spatially coherent over- 
densities, over a limited range of metallicities. To cre- 
ate the map we project the RGB stars onto a tangent 
plane, for which we chose a center at RA=0h 24 m 
and Dec = +44° which is near the center of the NW 
stream. This center is offset from M31's center, to make 
the co-ordinates reasonably rectangular in the region of 
the stream analysis. Without the RGB star selection no 
streams are readily visible in the halo. Figure [T] plots 
the sky distribution (in the tangent plane co-ordinates) 
of the [Fe/H] = [-0.6, -2.4] stars which we find below com- 
prise the stream. The same feature is clear in any map 
made with the low metallicity star sample. 

The low metallicity map, Figure [T] has a number of 
geometric sampling variations that must be controlled 
to make a reliable measurement of local stream density. 
The simplest problem is that the CCD array has gaps 
between the individual chips of ^80" in the horizontal 
direction and ~13" in the vertical direction. The images 
were dithered to cover the smaller gaps, but the resulting 
depth is shallower than the surrounding area and no sim- 
ple count correction procedure that did not increase the 
noise was found. These gaps are identified in the cam- 
era CCD co-ordinate system and masked out of the map. 
A more complicated geometric bias results from the cat- 
alog being built up from stacked local images in which 
the stellar photometry was done to create a local catalog. 
These catalogs were then matched in the areas of over- 
lap to eliminate duplicates and produce a global catalog. 
The images are vignetted towards the edge of the array 
so the image depth drops at the field edge, however the 
vignetting is partially compensated in the catalog match- 
ing at the image edges. That is, if the surface density is 
corrected with the inverse of the mean star density over 
the array then the map has an excess at the edges. It is 
straightforward to devise a correction which flattens the 
average surface density. The correction works well over 
a range of about 20% in surface density, but undesir- 
ably amplifies the noise if applied at lower completeness 
levels. We apply a cut to the image where the local sur- 
face density in CCD image co-ordinates drops below 20% 
of the average. The resulting masked out area is about 
15% of the image, a fraction which is not very sensitive 
to the precise value of the cut. The masked out regions 
are not included in the analysis and are readily visible in 
Figure]!] 

Fieldto field photometric depth variations will intro- 
duce artificial variations in the density of stars in the 
map. The depth of an image for a fixed exposure time 
depends on the sky transparency, brightness and image 
quality at the time of observation. Below the complete- 
ness limit there is substantial field to field variation that 
leads to a random checkerboard appearance, but the pat- 
tern essentially disappears above the completeness limit. 
Our measurement of the stream density subtracts the 
local mean background to reduce the problems of local 




Fig. 3. — A Mexican top-hat filtered image (inner radius 0.4° 
and outer 0.8°) that emphasizes the structure, but is not used for 
analysis. The image covers the same area as Fig.JT] with tick marks 
at 0.5°. The full width of the region used in the analysis, of which 
the outer half is used to define the local background, is shown. The 
gray scale extends from to 3 RGB stars per kpc 2 , greatly em- 
phasizing the contrast in the lumps. The elliptical angles, defined 
from the r(cj>) = a cost/) + bsin0 fit once rotated and centered, be- 
gin at <p=-75° (nearly horizontal line from the center) and proceed 
counter-clockwise in increments of 10° to </>=+70° (nearly vertical 
line). Note that these angles are not equal to the geometric angle. 
For our density variation measurements we use the well defined 
lower branch of the stream tj>= [—70, —10]. 

depth variations. 

A second, more difficult, form of depth variations is 
field to field variations in the zero points of the two filter 
bands. We can set a limit on these by using the double 
measurements of star brightness in the overlap regions. 
For bright stars where the photon noise is small the stan- 
dard deviation is 0.037 mag in g band and 0.052 mag in i 
band, with no discernible pattern over the field. The star 
catalog in our selected range of color and magnitude has 
nearly constant numbers with increasing depth near the 
completeness limit, so a 5% change in zero point leads 
to about a 5% change in numbers. The 5% should be 
considered an upper limit since the edges of the images 
are vignetted and have significant local calibration vari- 
ations. Across each 1° Megaprime field there are a myr- 
iad of camera co-ordi nate dependent photometric offsets 
( Regnault et al.|2009 ) which the current calibration does 
not fully take into account. These offsets lead to a change 
in zero point within each one degree pointing that is an 
additional ^2-3% variation. The masking procedure re- 
moves the worst of these variations, and the non-linear 
correction to flat star counts substantially reduces the 
variations within a field. What remains is below the shot 
noise of the star counts. 

2.1. The North- West Stream 

The North- West stream, readily visible in the star map 
of Figure[lj extends more than 100 kpc from M31 in pro- 
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TABLE 1 

[Fe/H] dependence of mean surface density 



[Fe/H] - 




[Fe/H] - 0.2 


stars kpc — 2 


-0.60 


0.05 


-0.80 


0.16 


-1.00 


0.41 


-1.20 


0.53 


-1.40 


0.27 


-1.60 


0.28 


-1.80 


0.06 


-2.00 


0.17 


-2.20 


0.08 


-2.40 


0.01 



jection. An ellipse is used to trace the stream, which cor- 
responds to finding the lowest order Fourier mode, that 
is, r((j>) = acos(0) + 6sin(0), rotated and centered for 
a best fit. Note that the angle </> is the angle once the 
ellipse is transformed into a sphere and is not the geomet- 
ric angle except along the major and minor axes. Higher 
order Fourier terms are not yet justified given the sig- 
nal to noise of the stream within these data. The ellipse 
parameters are varied, along with the stream width and 
the background region width to maximize the mean over- 
density within ±0.4° of the ellipse over the lower branch 
of the stream, taken as the angular range <\> = [—70, —10]. 
The stream region is shown in the Mexican top-hat fil- 
tered Figure[3] This filtered image, with an inner positive 
region of radius 0.4° and outer negative region of radius 
of 0.8°, with an area integral of zero, usefully illustrates 
the stream and noise properties near the filter length. 
The analysis is done on an unfiltcrcd high resolution ver- 
sion of Figure[T]with pixels of 18", ten times smaller than 
in Figures [T] and [3] 

We find that the best fit ellipse has a location of [- 
1.25°, +1.90°] in the tangent plane co-ordinates of Fig- 
ure ll|relative to our field center at [Oh 24m, +44°]. There 
is no particular reason for this fitting procedure to be 
constrained to have an ellipse centered on M31. The el- 
lipse has axes of 8.55° x 2.40° (or 117 x 33 kpc) at an 
angle of 42.25° W of N. The best fit occurs for an ellipse 
that is within the image data rather than for one that 
extends off the imaged area, even if the fitting procedure 
is started with a longer ellipse that extends out of the 
field. 

The straightforward consequences of projection are 
that the stream distance from M31 must be 117 kpc or 
larger, and, the ellipse must be at least as circular as 
we found above. We have assumed that all the stars in 
the stream are at 785 kpc, the adopted distance to M31. 
However, the distances of stars around the stream could 
vary roughly 100 kpc closer or further away, depending 
on how the stream is projected. Such distance variations 
will cause a distance modulus error of nearly ±0.3 mag, 
which means that individual stars will be identified with 
isochrones of the incorrect metallicity. This would be a 
significant problem if the stream were a narrow range in 
metallicity. However, our unweighted selection of stars 
over the very broad range of [Fe/H] = [-0.6, -2.4] corre- 
sponds to selecting a large rhombus in color-magnitude 
space that is 3.1 mag high in absolute magnitude. The 
upper limit to the distance differences of the stream rel- 
ative to the body of M31 would blur out our selection 
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Fig. 4. — The mean density transverse to the best elliptical fit 
to the stream over the (f> = [—70, —10] range. Errors are estimated 
from \Z~N in the bin. 

of stream stars about 10% or so, but it is only a small 
dilution of the stream density. 

2.2. Structure perpendicular to the Stream 

We measure the density perpendicular to the ellipse 
ridge line and display the result in Figure |4j The stream 
has a half width at half of the peak density of about 
2.5 kpc or a full width of 5 kpc. The stream width is 
dominated by the physical width of the stream, but also 
includes a significant component from the stream loca- 
tion deviating slightly from our assumed elliptical shape. 
The stream subtends an angle ranging from about 0.1 
radian to 0.05 radian, as seen from the center of M31, 
although the stream is physically wider than the cool, 
but relatively close, streams in our own galaxy. 

The location of the centerline of the stream, calculated 
as the density weighted mean distance from the best fit 
ellipse, is shown in Figure [5j The error flags are calcu- 
lated from the statistics of the star counts and confirms 
the impression of Figure [3] that the stream center varies, 
sometimes discontinuously, from the best fit ellipse. The 
variations are largest on the upper branch of the stream, 
4> > 0, where it has large gaps. On the more continuous 
lower branch there is a statistically significant offset from 
the mean centerline at <fr — —30°. 

2.3. Metallicity and Luminosity Distribution of the 

Stream 

The density is measured within a constant width re- 
gion around the ellipse, ±0.4° (~ 5 kpc). We subtract 
the foreground/background measured at the same angu- 
lar bin within the two adjacent 0.2° strips on either side 
of the stream. The range of metallicity in the stream 
is shown in Table [T in bins of width 0.2 in [Fe/H]. Ta- 
ble [T] gives the top of the metallicity bin in column 1, the 
mean number of stars per kpc 2 in column 2. This table 
shows that the stream has its maximum density in the 
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Fig. 5. — The location of the stream center relative to the best 
fit ellipse centerline, for 2° and 4° binning. The variations are 
larger than the expected statistical variations, particularly for the 
"tilted" segments at positive angles. 



Fig. 6. — The mean density around the stream in bins of 2° of 
elliptical angle. The errors are estimated from the star counts in 
the bins and the surrounding background region. The \ 2 statistics 
is calculated only for the <j> = [— 70, —10] range. 



[Fe/H]~ — 1.2 bin. For our measurement of the stream 
density we will use the entire [Fe/H] = [-0.6,-2.4] range, 
uniformly weighted. In principle somewhat better sig- 
nal to noise would be obtained if we weighted the stars 
with the metallicity distribution function. However, the 
flat distribution gives the most straightforward stream 
density measurement. 

The stream is broadly distributed in brightness with 
the mean number of stars [21-22, 22-23, 23-24, 2x(24- 
24.5)] i-band mag in the [Fe/H] = [-0.6, -2.4] range be- 
ing [0.62, 0.78, 0.55, 0.47] per magnitude per kpc 2 over 
the -70 to -10 angle range. The luminosity function of 
the stream is essentially flat over the absolute magnitude 
range Mj < +0 as is usually the case for low metallic- 
ity RGB luminosity functions. Physical properties of the 
stream are derived below. 

3. DENSITY STRUCTURE OF THE STREAM 

Whether the NW stream is smooth or lumpy is the 
observational question of primary interest in this paper. 
This measurement needs to be handled carefully since 
the stream over-density is only about 15% of the total 
nominal RGB star density. The stream density will vary 
simply due to the \/~N fluctuations in the few thousand 
stars in each angular bin. The density measurement is 
done within the region shown in Figure [3j The stream 
is defined as being a fixed width of the centerline and 
the local background is measured within the cutout re- 
gion beyond the stream. The analysis region of both the 
stream and local background is shown in Figure [3| The 
ne wly id entified dwarf galaxy AndXXVII ( jRicharctson et 



al 



2011), just above the upper branch of the stream in 
eRj is masked out to exclude it from the local back- 
dcorrection. 

The density along the stream, within a half-width of 
0.4° on the sky subtracting the local background in the 




Fig. 7. — Same as Fig. [3] showing that the stream is not evident 
in the [Fe/H] = [-0.5, 0]. The contrast scale is the same as in Fig. [3] 

adjacent 0.2° regions, is shown in Figure [6] at intervals 
of 2° in the ellipse angle, (j>. The errors are computed 
from the y/~N of the star-counts in the stream and in the 
background region, added in quadrature. 

3.1. Statistical Significance of Density Variations 

The x 2 test estimates the probability that a set of data 
is consistent with a model, which in our case is a constant 
density. We need to be aware that sampling and calibra- 
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TABLE 2 
Probability of a Smooth Stream 
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Fig. 8. — The angular dependence of the stream density at high 
metallicity, showing that the density variations due to background 
are not significant. The error flags at each point are about 30% 
smaller than for the stream because there are fewer stars. 



tion variations could artificially create density variations. 
Masking takes care of the CCD gaps and the low expo- 
sure regions of the images. The field flattening largely 
corrects the counts over the rest of the image. Any resid- 
ual photometric fluctuations will be measured in the con- 
trol fields to estimate their size. 

The x 2 statistic applied to the </> = [—70,-10] range 
indicates that the chance that the stream has a uniform 
density is less than about 10 -6 . A binning angle of 2° 
corresponds to roughly 2-4 kpc, depending on location 
around the ellipse. We have measured the stream den- 
sity variations in bins ranging from 1° to 8° finding that 
the variations are present at all these scales. Since there 
are easily visible clumps and gaps for <j) > — 10 including 
the upper branch of the stream in the % 2 analysis would 
significantly increase % 2 , which of course could well be a 
real description of the larger stream. We prefer to present 
the most conservative statistic on the visibly smoothest 
part of the stream. We emphasize that even after re- 
stricting our statistical analysis to the lower branch, the 
segment is the longest stream segment available for study 
so far. Moreover its orbit appears to keep it well away 
from the potential disturbances of the structure in the 
disk of M31. 

We have varied virtually all stream and image param- 
eters to check the result that the density variations are 
truly those of the star stream. Reasonable changes in the 
path of the ellipse, its width and the local background 
make little difference, other than slightly reducing the 
mean density of the stream. The photometric zero point 
variations that remain in this version of the data are 
readily estimated. The mean background is approxi- 
mately 12 stars kpc -2 and the photometric variations 
are expected to produce 5% or smaller fluctuations co- 
herent over scales of 1 degree in the image, or a standard 
deviation of about 0.6 stars kpc -2 . The measured den- 



l-P(xV) 



NW Stream 

1 2.18 59 4 X 10" 

2 3.01 28 2 x 10" 
4 3.58 14 6 x 10" 



Null test with [Fe/H]>-0.5 

1 1.17 59 0.17 

2 1.16 28 0.25 
4 0.65 14 0.85 

2° X 4° offset stream 

1 1.96 29 0.002 

2 1.57 13 0.090 
4 2.38 6 0.027 



sity fluctuations have a standard deviation of ~2.0 stars 
kpc -2 , at least three times the size of the variations that 
the photometric zero points are expected to cause. 

A straightforward control sample is to measure the 
fluctuations at the same location in higher metallicity 
stars, [Fe/H]>-0.5, in which the stream is not detected. 
Figure [7] shows the Mexican top-hat filtered image in 
the hign metallicity stars, demonstrating that the back- 
ground fluctuations have a significantly different pattern 
but have a similar amplitude. The density around the 
stream is displayed in Figure [§] The x 2 test for the sig- 
nificance of the density variations is reported in Table [2} 
No significant variations are found. 

A second control sample is to simply offset the mea- 
surement region away from the stream. We add 2° to 
the minor axis and 4° to the major axis and repeat the 
measurements, adjusting the (f> range to approximately 
the same distance. The statistics are reported in Table [2] 
finding that the probability that the variations are due to 
chance to be about 5%. Moving the region around find 
finds that the \ 2 values fluctuate a fair bit depending 
on position although nothing as significant as the NW 
stream itself ever emerges. These fluctuations are likely 
the result of shorter and lower significance star streams 
that are real entities. 

3.2. Upper Branch Stream Substructure 

Whether the upper branch of the stream in Fig. [I] or [3] 
is the extension of the lower branch of stream is unclear 
in these data. The upper branch certainly fits within the 
same ellipse and has a shared metallicity distribution at 
the precision we can measure it. However, the upper 
branch is less a stream than three dominant segments 
separated by gaps of comparable size. It is of great in- 
terest to note that the segments are tilted with respect 
to the general path of the stream. Simulations show that 
tilted segments are some of the dominant structu res fre- 
quently seen in chopped streams, as first noted in |Yoon,| 
Johnston & Hogg (2010), and as we show below. 



3.3. Dynamics and Alternate Clumping Mechanisms 

The total luminosity of the lower branch of the stream 
is found by summing the over-densities in one magnitude 
bins within a range of [-70,-10] ecliptic angle, to find a 
stream segment luminosity of 4.3 x 10 5 L Q to Mj ~ mag. 
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The luminosity is the sum over the RGB stars alone. To 
estimate the total luminosity we use the cumulative stel- 
lar luminosity function of a stellar population of similar 
age and metallicity, for which we chose the globular clus- 
ter M12. M12 is an intermediate metallicit y cluster with 



a wel l studied stellar luminosity function ( |Hargis et al. 
2004|. Using their distance modulus of 14.0 mag and 



their I band luminosity function we find about 58% of 
the light above Mi = mag. Therefore the corrected 
total luminosity is 7.4 x 1O 5 L . 

The mean surface brightness is 1060 L kpc -2 . If we 
assume the stream is a uniform density cylinder of 5 
kpc radius and a stellar mass-to-light ratio of, say, 3 
in the i-band, then the volume stellar mass density is 
500 M kpc -3 . This can be compared to the mass den- 
sity in the dark matter halo at say 90 kpc, which for 
our adopted halo parameters is 1.2 x 10 5 M kpc -3 , that 
is, about a factor of about 240 higher than the stellar 
stream mass density. At 30 kpc the ratio is 3600 times 
more background halo dark matter than stream stellar 
mass. If dark matter from the progenitor is mixed into 
the stream then the mass density of the stream would 
be proportionally increased, lowering the local mass ra- 
tios but the stream is unlikely to be significantly self- 
g ravitating anywhere. 



Quillen & Comparetta (2010) have raised the possibil- 
ity that Jeans instabilities can be the source of of density 
variations in the stream. The Jeans length in the stream, 
A j = \J~ko 2 /Gp evaluated with the stellar mass density 
and assuming an internal velocity dispersion of 10 km s -1 
is about 400 kpc. On the other hand if the velocity dis- 
persion were as low as 1 kms -1 then the formal Jeans 
length would be 40 kpc, comparable to the longest scale 
of density variations in the stream. Very low velocity 
dispersions are expected for slowly tidally stripped glob- 
ular clusters. The presence of several globular clu sters 
along the course of the stream ( |Mackey et al.||2010 ) and 
the possibility that AndXXVII ( |Richardson et al.|[20TT | 
is the source of the stream suggests that the internal ve- 
locity dispersion is at the level that a dwarf galaxy would 
create, ~ 3kms -1 or more. The resulting Jeans length 
is then roughly 120 kpc, the size of the stream. It is 
important to note that the stability analysis is far from 
complete. First, the disrupted object may have its own 
dark matter which has identical kinematics to the visible 
stars and would increase the mass density. Second, the 
stability analysis needs to account for the very strong 
tidal field of the M31 dark halo. Strong tides generally 
act to suppress much instability, particularly in largely 
radially oriented streams, which is likely the case for the 
NW stream. 

Another mechanism that generates clumps in tidal 
streams is simply the pil eup of stars at low ve locity points 
of their epicyclic orbits ( Kiipper et al.|200"8 ). The lumps 
appear at 127r times the tidal radius. Since no progen- 
itor object for the stream has been confidently identi- 
fied, the tidal radius is very uncertain. However if we 
take a low luminosity dwarf galaxy of stellar plus a rel- 
atively low dark mass for a total progenitor dwarf mass 
of Md — 10 s M , then at a galactic radius of r g = 100 

kpc, the tidal radius, r s (Mrf/(3M s )) 1 / 3 , is approximately 
3.2 kpc. The epicycle pileup separation will then be at 



a separation of 120 kpc. That scale compares to the 
size of the stream. The stellar mass in the stream alone, 
2 x 10 6 M for a fairly minimal assumed mass-to-light 
ratio of 3, gives about 40 kpc epicyclic clumping scale, 
which is at the upper end of range of what we see, but 
still leaves the smaller scale variations to be explained. 
We conclude that epicyclic pile-up is not a significant 
clumping factor for the NW stream. 

If the progenitor system had large internal density vari- 
ations in a highly order velocity field they could be fed 
out into the stream to create lumps. However, this idea 
runs into trouble with our astrophysical knowledge of low 
mass galaxies. Dwarf spheroidal galaxies are comprised 
largely of stars on randomly oriented orbits and have es- 
sentially no substructure, which when tidally disrupted 
would lead to a smooth stream. Dwarf irregular galaxies 
do have substantial substructure in their gas and young 
stellar populations. The visible luminosity (corrected for 
light below our observational limit) in the lower branch of 
the stream is only 7 x 10 5 L . Dwarfs with such low lumi- 
nosities usually have comparable rotational and random 
velocities which would quickly blur out stellar structures 
deposited into a stream. Given that |Mackey et aL (2010) 
associate 3 or 4 globular clusters with the 1NW stream, 
which would lead to a unusually high ratio of globular 
cluster luminosity to progenitor luminosity, it it possible 
that the NW stream is only part of a larger stream that 
has yet to be clearly identified. The structures within 
star forming dwarf galaxies are often low mass stellar as- 
sociations and clusters that fairly quickly dissolve as a 
result of stellar mass loss and their own internal dynam- 
ics. Spiral patterns are at best very chaotic in low mass 
dwarf irregulars. Our assessment is that although care- 
fully sub-structured young progenitor disk system could 
produce some of the stream lumpiness, appropriate as- 
trophysical systems do not exist. 



4. COOL STREAM SIMULATIONS 

Star streams normally originate from the tidal disso- 
lution of cither dwarf galaxies or globular clusters. In 
both cases tidal fields pull off outer stars or sometimes 
dissolve the entire system and distribute the stars in 
nearby orbits around the galaxy. The effect of adding 
randomly orbiting dark matter sub-halos on hi ghly ide- 
alized star streams has been studied prev iously fllbata et| 
aLl|2002l [STegal-Gaskins fc VaUuri| |2508l |CarlbergP(jU9t 
Yoon, Johnston & Hogg 2010), finding that the streams 



were folded, chopped and heated to about 30kms over 
a Hubble time. For the examination of the expected de- 
gree of lumpiness of star streams we need a more realistic 
model stream. T he association of globular clusters with 
the NW stream ( |Huxor et aL|2008l packey et al.||2010 1 
indicates that a dwarf galaxy, likely AndXXVII, is the 
source of the stream. 

The study of the effects of dark matter sub-halos on 
the stream is inherently statistical since we are unsure 
to what degree the visible stars and gas give a complete 
census of the sub-halos, as highlighted in the celebrated 
missi ng satellite problem ( Klypin et al.|1999||Moore et al. 
1999 1. A detailed match to the orbit ot the N W stream is 



not required for this very basic assessment of the action of 
sub-halos on a stream. The details of the stream creation 
process have relatively little to do with the subsequent 
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sub-halo interactions with the stream. Our test particle 
stream originates from particles orbiting in a constant 
mass Plummer sphere to which particles in the core re- 
gion are gradually given velocities that boost them out 
into the region where tides can carry them away. Future 
studies will report a full suite of self-gravitating simula- 
tions. 

The simulation uses a galactic potential with a spher- 
ical NFW halo pote ntial scaled to the re sults of the 
Aquarius simulation ( |Springel et al. 2008). The rota- 
tion curve of M31 is somewhat higher than the Milky- 
Way, but beyond 30 kpc mass modeling finds that th e 



two galaxies are very similar (Carignan et al. 2006) 



For simplicity and comparability we therefore ado pt~the 
same disk-b ulge potential for the disk and bulge as 



ston 



John- 



( 1998|) which has a M iyamoto-Nagai disk potential 

J iJe_ "NT :llinTch — j „ 1 — 



I Miyamotd& iNagai||l975l) and a nuclear bulge modeled 



with a Plummer sphere. The sub-halos are modeled as a 
collection of Plummer spheres following the distribution 
in masses, o rbits and internal structure of the Springcl 
et al. (2008) simulations. The star stream particles are 
evolved in the combined potentials of the background 
galaxy and orbiting sub-halos. We use time steps of 
0.14 Myr, or about 10 5 in a Hubble time, to ensure that 
the quickly moving particles are integrated properly over 
the small scale potentials of the sub-halos. 

4.1. Density Variations in the Simulations 

The simulations start with the Plummer sphere 10 6 M Q 
progenitor located at x = 100 kpc from the center in the 
plane of the disk, with y = z = 0. The initial velocity 
is chosen to be entirely tangential perpendicular to the 
plane of the disk with an angular momentum of about 
50% of the circular value at that radius. In the absence 
of sub-halos there are no density variations in the simu- 
lated star streams. M31 contains nearly 30 known dwarf 
galaxies. To what degree could these dwarfs, along with 
their dark-halos, cause the density variations we are mea- 
suring? For instance, what is the effect of the heaviest 
5 sub-halos alone? Figure [9] shows the x — y plane pro- 
jection of the stream in the presence of 0, 5, 100, and 
1000 sub-halos as drawn from the M~ 19 mass distribu- 
tion. Including the 5 most massive sub-halos produces 
only large scale distortions of the stream relative to the 
smooth halo. With 100 sub-halos we see significant small 
scale variations and going to 1000 sub-halos adds yet 
more small scale structure. Figure [l0| shows the density 
variations in angular bins assuming that the stream has 
a constant width, which is close to being equivalent to 
the M31 measurements shown in Fig. [6l In this partic- 
ular simulation the progenitor remains bound, so there 
is a density peak near <fi = 0. Away from the progeni- 
tor there is a clear trend of increasing disturbance in the 
density variations as the number of sub-halos increases. 
Neither a smooth halo nor one with only the five most 
massive sub-halos, which roughly mimics the most mas- 
sive dwarf galaxies, produce anywhere near the amount 
of substructure we see. One hundred sub-halos do a rea- 
sonable job in the outer parts of the tail but do not do 
much damage to the inner part of the tail. Once we get to 
one thousand sub-halos we star seeing gaps and clumps 
essentially everywhere. 

Deviations from a mean centerline are also a signal of 
sub-halo interactions. Again the hundred halo simula- 
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Fig. 9. — An example of stream disturbances under the action 
of 0, 5, 100 and 1000 dark matter sub-halos at an age of 11.12 
Gyr. Note the increasing small scale variation as the number of 
sub-halos increases. 
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Fig. 10. — The density profiles in 2° angular bins around the 
streams of Figure[9] Away from the progenitor at = (which has 
remained bound in this case) the density variations approximately 
mimic the character of the M31 stream. Note that the stream is 
double valued around 90°. 



tion produces a few angled clumps but they are much 
more clearly present with the thousand halo simulation. 
We will report on the statistics of a much larger set of 
simulations elsewhere. 

In this illustrative study we find that of order of one 
thousand sub-halos are required to produce all the small 
scale structure on the scales see n in the M31 NW stream 
from a n initially smooth stream. Yoon, Johnston &: Hogg| 
(20101 report a more complete dynamical analysis fo- 
cused on streams at smaller radii, but the dynamical out- 
comes they describe are quite generally applicable. They 
also firmly conclude that a large population of dark mat- 
ter sub-halos will induce gaps and clumps in streams. A 
larger statistical study will be reported elsewhere. 

5. DENSITY STRUCTURE OF OTHER STREAMS 

A spe ctacular star stream in ou r own galaxy is the Pal 5 
stream ( Odenkirchen et al7||2001 ) which i s kinematically 
cold and narrow ( Oden kirchen et al.|2 009 ) . At a distance 
of 23.2 kpc the visible part is only 4 kpc long. It is very 
narrow, about 0.4 kpc, roughly 1/5 of the width of the 
M31 NW stream segments. The source of the stream 
is the outer halo globular cluster Pal 5 with the stream 
having the highest surface density a t the location of its 
rogen itor. The density analysis of Odenkir chen et al.| 
2003[ ) finds that there are density variations of the Pal 
stream above the shot noise in the northern part of the 
stream, but not the southern part. The significance of 
the variation is not given, but there are 9 of 25 points 
more than la from the mean trend over about 1 kpc 
of length. This certainly constitutes the detection of a 
significant lump in the currently available 4 kpc length 
of the Pal 5 stream, but on a smaller scale. 



Johnston, Spergel & Haydn ( 2002 ) were inspired by the 
narrow velocity and geometric distribution of the 75 car- 



bon stars of the Sagittarius stream ( Ibata et al. 12001' ) to 
consider the lumpiness of the Milky Way halo. They de- 
velop a statistic based on the fifth through tenth Fourier 
terms to compare smooth halo models to models with up 
to 256 sub-halos, being careful to consider models that 
allow for the main dwarf galaxies. When applied to the 
available Sagittarius data they find that "the degree of 
scattering is entirely consistent with debris perturbed by 
the LMC alone." If applied to the 100 kpc NW M31 
stream, the tenth Fourier term would be a 0=36° varia- 
tion, which is roughly comparable to a 16° binning where 
there still is real structure relative to a constant, but only 
o n yet larger scale s. 



Majewski et al. 



( |2004[ ) considered kinematic data for 
the Sagittarius stream m our galaxy and concluded that 
the halo could not be very lumpy, however we note 
that our simulations show that sub-halos tend to "chop" 
streams leaving locally cold remnants with pieces at the 
same radius offset in velocity, wh ich may be compati- 
ble with the Majewski et al. (2004) Figure 2. Kinematic 
tests will become very powerful in testing for substruc- 
ture once large velocity samples are available. 

6. DISCUSSION AND CONCLUSIONS 

The 100 kpc NW stream in M31 is a coherent geo- 
metric structure whose orbit is well clear of the body of 
M31 making it a near ideal testing ground for the pres- 
ence or absence of thousands of dark sub-halos predicted 
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in LCDM n-body simulations. Although the stream ap- 
pears to be nearly a half ellipse over a common range of 
metallicities the upper branch of the stream is less well 
defined and has a number of clearly visible gaps. The 
lower branch is nearly complete and provides a much 
more conservative test for the presence of sub-halos. The 
main result of this paper is that the stream has highly 
significant density variations on virtually all scales from 2 
kpc, up to about 20 kpc. The variation of density around 
the mean has a very low probability of being a chance sta- 
tistical fluctuation, less than 10~ 5 . We have been careful 
to measure the lumpiness relative to an averaged local 
background and masked out gaps and regions where the 
photometric uncertainties add significantly to the varia- 
tions. As a control sample we take a higher metallicity 
set of stars at exactly the same location as the stream, 
which finds no significant density variations. Relative to 
other known cool star streams, the M31 stream stands 
apart for its length, distance from the disturbing effects 
of the M31 disk, a variety of scales of substructure and 
its high statistical significance. 

It is interesting to note that M31 appears to have only 
the one well defined long stellar stream at large radius. 
The NW stream has quite a low total luminosity so it 
should not particularly stand out relative to other stel- 
lar streams that were created at large radius over the 
buildup of the M31 dark halo. Several other stream frag- 
ments are visible in the full field map Richar dson et ah] 



( 2011 1 but none as long or coherent as the NW stream. If 



one assumed that other long streams likely formed over 
the lifetime of the halo then the absence of others at 
the present time indirect ly suggests tha t they have been 
broken up by sub-halos ( |Carlberg||2009[ ) . 

The measured NW stream lumpmess essentially rules 
out the possibility that it is a low mass star stream or- 
biting in a smooth galactic potential of a disk plus bulge 
plus dark halo. Conversely, the NW stream density vari- 
ations are compatible with the level of density changes 
that a large population of dark matter sub-halos induce. 
The details of the density variation are sensitive to the 
statistical distribution of the sub-structure but are un- 
likely to be specifically modeled for a single stream. As 
more streams are observed statistical modeling of the 
degree of substructure present should become possible. 
The present study is statistically consistent with a highly 
sub-structured dark halo, although we caution that the 
density data alone for a single stream is not a conclusive 
proof that halos are as sub-structured as LCDM simula- 
tions predict. Overall, the density variations of the NW 
stream are strong circumstantial evidence that the pre- 
dicted thousands of dark matter sub-halos are present in 
M31's dark halo. 
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